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A NEW MECHANISM FOR CORE-COLLAPSE SUPERNOVA EXPLOSIONS 

A. BuRR0ws\ E. LivNE^, L. Dessart\ CD. Ott^, J. Murphy^ 
Accented to Ap.J. 

ABSTRACT 

In this paper, we present a new mechanism for core-collapse supernova explosions that relies upon 
acoustic power generated in the inner core as the driver. In our simulation using an 11-Mq progenitor, 
an advective-acoustic oscillation a la Foglizzo with a period of ^25—30 milliseconds (ms) arises ^^200 
ms after bounce. Its growth saturates due to the generation of secondary shocks, and kinks in the 
resulting shock structure funnel and regulate subsequent accretion onto the inner core. However, this 
instability is not the primary agent of explosion. Rather, it is the acoustic power generated early on 
in the inner turbulent region stirred by the accretion plumes, and most importantly, but later on, by 
the excitation and sonic damping of core g-mode oscillations. An £ — 1 mode with a period of ^3 ms 
grows at late times to be prominent around ~500 ms after bounce. The accreting protoneutron star 
is a self-excited oscillator, "tuned" to the most easily excited core g-mode. The associated acoustic 
power seen in our 11-Mq simulation is sufficient to drive the explosion >550 milliseconds after bounce. 
The angular distribution of the emitted sound is fundamentally aspherical. The sound pulses radiated 
from the core steepen into shock waves that merge as they propagate into the outer mantle and deposit 
their energy and momentum with high efficiency. The ultimate source of the acoustic power is the 
gravitational energy of infall and the core oscillation acts like a transducer to convert this accretion 
energy into sound. An advantage of the acoustic mechanism is that acoustic power does not abate 
until accretion subsides, so that it is available as long as it may be needed to explode the star. This 
suggests a natural means by which the supernova is self-regulating. 

Subject headings: supernovae, neutrinos, multi-dimensional radiation hydrodynamics, stellar pulsa- 
tions 



1. INTRODUCTION 

The essence of the mechanism of core-collapse super- 
novae must be the conversion of a fraction of the reservoir 
of gravitational energy into the kinetic and internal en- 
ergy of the exploding mantle of the Chandrasekhar core 
whose instability inaugurates core collapse. A protoneu- 
tron star (PNS) is left which evolves into a cold neutron 
star most of the time, or, rarely, is the intermediate state 
along the way to the formation of a stellar-mass black 
hole through the general-relativistic instability Due 
predominantly to prodigious neutrino radiation (losses) 
over hundreds of milliseconds to seconds, the PNS be- 
comes more and more bound, on its way to achieving a 
total binding energy near ~3xl0^^ ergs, a mass-energy 
equivalent of ~0.15 Mq. This is far in excess of the 
^10^"'^ ergs of the supernova explosion, so it has long 
been thought that 1) the mechanism of explosion is very 
inefficient, and 2) there is so much energy in the neutrino 
emissions that if only a small fraction were captured in 
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the mantle it could be driven to explosion. The latter 
notion is the essence of the "neutrino-driven" mecha- 
nism, first adumbrated in its direct form by Colgate & 
White (1966) and Arnett (1966) and in its delayed form 
by Wilson (1985) and Bethe & Wilson (1985). However, 
using the best input physics and numerics, it has been 
shown that in spherical symmetry both the prompt and 
the delayed-revitalization neutrino heating mechanisms 
fail (Rampp & Janka 2000; Licbcndorfer et al. 2001; 
Thompson, Burrows, & Pinto 2003). Moreover, the sug- 
gestion by Wilson (1985), Wilson & Mayle (1988,1993), 
and Mayle & Wilson (1988) that "neutron-finger" insta- 
bilities interior to or near the neutrinospheres can advect 
heat, and thereby boost the driving neutrino luminosi- 
ties sufficiently to lead to pseudo-spherical explosions, 
has been dealt a severe blow by the more detailed analy- 
ses of Bruenn & Dineva (1996), Bruenn, Raley, & Mezza- 
cappa (2004), and Dessart et al. (2005). Finally, through 
two decades of simulations, the idea that on dynamical 
times the bounce itself could, by direct piston action on 
the outer mantle, launch a successful shock wave that 
does not stall into accretion has been refuted. Due to 
significant shock breakout neutrino losses and nuclear 
photodissociation, the shock must stall. 

Out of the ashes of the one-dimensional (ID), spher- 
ical models was born the idea that multi-dimensional 
neutrino-driven convection could increase the efficiency 
of neutrino heating behind the stalled shock and lead to 
a delayed explosion (Hcrant et al. 1994; Burrows, Hayes, 
& Fryxell 1995; Janka & Miiller 1996; Fryer & Warren 
2002,2004). Such neutrino-driven convection arises nat- 
urally and does increase the radius of the stalled shock 
and the size of the so-called "gain region" external to 
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the neutrinospheres where there is net neutrino energy 
deposition. This is the current paradigm and models that 
did not explode in ID have indeed exploded in 2- and 3- 
D. However, all these previous successful explosions were 
obtained in the context of simplified or gray treatments 
of neutrino transport and it has long been suspected 
that more sophisticated, multi-group, multi-angle treat- 
ments of neutrinos in the multi-dimensional spatial con- 
text were needed. In fact, though implementing 2D neu- 
trino transport in an approximate way, Rampp & Janka 
(2002) and Buras et al. (2003,2005) use 2D hydrodynam- 
ics and one of the best available multi-group co-moving- 
frame Boltzmann solvers and do not obtain explosions'". 
Despite incorporating the best input physics, and coming 
tantalizingly close to explosion, their simulations lead to 
fizzles, calling into question the neutrino-driven mecha- 
nism itself. Within 200-300 milliseconds of bounce, after 
growing out to 150-200 kilometers, their shocks recede 
and are declared to be unsuccessful. 

If the energy transfer from the core to the mantle nec- 
essary to explode the star is by neither direct hydrody- 
namics nor neutrino heating, what is left? How do core- 
collapse supernovae explode? In this paper we propose a 
new alternative, the generation in the core and the prop- 
agation into the mantle of strong sound waves. Acous- 
tic power is, potentially, an efficient means to trans- 
port energy and momentum into the outer mantle to 
drive the supernova explosion. Unlike neutrinos, sound 
is almost 100% absorbed in the matter, though some 
of the acoustic energy is reradiated by neutrinos. As 
sound pulses propagate outward down the density gra- 
dient they steepen into multiple shock waves that catch 
up to one another and merge; a shock wave is almost a 
perfect black-body absorber of sound. If sufficient sound 
is generated in the core, it would be a natural vehicle 
for the gravitational energy of infall to be transferred to 
the outer mantle and could be the key missing ingredi- 
ent in the core-collapse explosion mechanism. Further- 
more, periodic shocking due to multiple sound pulses can 
lead naturally to entropies in the debris of hundreds of 
units, just what is required for r-process nucleosynthesis 
(Woosley & Hoffman 1992; Woosley et al. 1994; Hoffman 
et al. 1996). 

We believe we have identified a vigorous source for 
the necessary acoustic power: the excitation and oscil- 
lation of core pulsation modes in the deep interior of 
the PNS. Using the 2D radiation/hydrodynamic code 
VULCAN/2D we have discovered that turbulence and 
anisotropic accretion in the inner 40—100 kilometers can 
excite and maintain vigorous core g-mode oscillations 
which decay by the radiation of sound. The inner core 
acts as a transducer for the conversion of accretion grav- 
itational energy into acoustic power. The associated 
acoustic power seen in our simulations is sufficient to 
drive the explosion >550 milliseconds after bounce. 

In we summarize the input physics, computational 
capabilities, and numerical approaches incorporated into 
the current version of our radiation/hydrodynamics code 
VULCAN/2D. Then, in i^lwe describe our baseline sim- 
ulation results for the 11 Mq progenitor star on which 
we focus. This section defines the various hydrodynamic 

^ unless they drop their many velocity-dependent terms, some- 
thing they do not advise doing 



stages and events and explores the important, but sec- 
ondary, role of the advective/acoustic or "SASI" insta- 
bility highlighted by Foglizzo and collaborators (Foglizzo 
2000,2001; Foglizzo & Tagger (2002); Foglizzo, Galletti, 
& Ruffert 2005) and by Blondin, Mezzacappa, & De- 
Marino (2003). We also introduce the key role played 
at late times by acoustic power and the emergence of 
core oscillation modes and discuss some aspects of the 
components of the mechanism we envision. Then, in SQI 
we summarize the general advantages of acoustic driv- 
ing, and follow this in f|5lwith an investigation into the 
excitation of core oscillation modes, the associated gen- 
eral physics, and their acoustic damping. After this, we 
summarize in fjSlthe emerging role of neutrinos as impor- 
tant, but secondary, players in the explosion mechanism 
we are presenting. Finally, in fJZ|we reiterate the central 
features of our new mechanism, explain why it was not 
discovered earlier, and speculate on what now needs to be 
done both to test our hypothesis and to further advance 
our understanding of this central problem in theoretical 
astrophysics. 

2. COMPUTATIONAL APPROACH AND INPUT PHYSICS 

The code VULCAN/2D uses the hydrodynamic ap- 
proach described in Livne (1993), with the transport 
methods discussed in Livne et al. (2004) and Walder 
et al. (2005). It is a Newtonian, 2D, multi-group, multi- 
angle radiation/hydrodynamics code with an Arbitrary- 
Lagrangian-Eulerian (ALE) structure (with remap) and 
a Multi-Group, Flux-Limited Diffusion (MGFLD) vari- 
ant. Velocity terms in the transport sector, such as 
Doppler shifts, are not included in the code, though ad- 
vection is. Note that the velocity terms in Eulerian trans- 
port are different from the coresponding terms in the co- 
moving frame and that general statements about their 
relative importance are very frame-dependent^. We par- 
allelize in energy groups using MPI and do not include 
redistribution by neutrino-electron scattering in the 2D 
code. Such redistribution and scattering are of mod- 
est import on infall, affecting the trapped electron frac- 
tion (Ye) and entropy (5) by ^10%, but are otherwise 
quite subdominate. At a neutrino energy of 10 MeV, the 
neutrino-electron scattering cross section is ~100 times 
smaller than the dominant cross sections off nucleons. 
Thompson, Burrows, & Pinto (2003) have devised and 
implemented an explicit scheme to handle redistribution 
that adds only '-^10 — 15% to the computational time 
and this approach will be incorporated in future versions 
of VULCAN/2D. Note that the attempt to handle the 
full energy/angle redistribution problem implicitly has 
resulted in codes that are thereby slower by many fac- 
tors (not percent), severely inhibiting their use for explo- 
rations in supernova theory. 

In this paper, due to its relative speed, we use the 
MGFLD implementation of VULCAN/2D. The flux lim- 
iter is a vector version of the one found in Bruenn (1985). 
The code can handle rotation, though the results re- 
ported here are for a non-rotating progenitor. In 2D, 
the calculations are axially/azimuthally symmetric, and 
we use cylindrical coordinates (r and z), but the grid 
points themselves can be placed at arbitrary positions. 
This allows us to employ a Cartesian grid at the center 

^ In particular, even the sign of the effect can change. 
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(inner ^20 kilometers) and transition to a spherical grid 
further out. The grid resolution is essentially uniform ev- 
erywhere within ~20km. A version of this grid structure 
is plotted in Ott et al. (2004). The Cartesian format in 
the interior allows us to avoid the severe Courant prob- 
lems encountered in 2D by all other groups employing 
grid-based codes due to the inner angular Courant limit 
and, thereby, to perform the calculations in full 2D all 
the way to the center. In many simulations to date, the 
inner core was calculated in ID and grafted onto an outer 
region that was handled in 2D (e.g.. Burrows, Hayes, & 
Fryxell 1995; Swesty & Myra 2005ab), or was excised 
completely (e.g., Janka & Miiller 1996; Blondin, Mezza- 
cappa, & DeMarino 2003; Scheck et al. 2004). In others, 
the 2D simulation was done on a ~90° wedge, with the 
inner 2 or 25 km done in ID (Buras et al. 2005). The 
gray SPH simulations of Herant et al. (1994) and Fryer 
& Warren (2002,2004) are an exception. Originally, a 
major motivation for this global 2D feature was the self- 
consistent investigation of core translational motion and 
neutron star kicks. However, as we will see, freeing the 
core has other advantages. 

To determine the gravitational potential field we use 
a finite-difference approach to solve the Poisson equa- 
tion. In addition, we incorporate the gravitational force 
along the symmetry axis in an automatically momentum- 
conserving fashion by writing it in divergence form. 
VULCAN/2D has an option to move the grid after 
bounce to maintain the best zoning under the core, 
whether it moves or not, while at the same time tracking 
this core motion. This feature ensures that the highest 
resolution is placed under most of the mass. Outside of 
the Cartesian mesh, the baseline calculation employs 121 
angular zones equally spaced over 180°, and logarithmi- 
cally allocates 162 radial shells between ^20 km and the 
outer radius at 3800 km^. 

We employ the equation of state (EOS) of Shen et 
al. (1998), since it correctly incorporates alpha particles 
and is easily extended to lower densities. The neutrino- 
matter interaction physics is taken from Thompson, Bur- 
rows, & Pinto (2003) and Burrows & Thompson (2004). 
The tables generated in T/p/ Ye /neutrino-species space 
incorporate all relevant scattering, absorption, and emis- 
sion processes. We follow separately the electron neu- 
trino (z/e) and anti-electron neutrino (j^e), but for com- 
putational efficiency we lump the four remaining known 
neutrinos into "i/^" bins in the standard fashion. Our 
baseline model has 16 energy groups for each species, 
distributed logarithmically from 1 to 200 MeV. More en- 
ergy groups would be better. Due to extreme matter- 
suppression effects, we have not felt it necessary to in- 
corporate the effects of neutrino oscillations (Strack & 
Burrows 2005). 

The instabilities that develop in the early stages of 
the post-bounce phase are seeded by the slight pertur- 
bations introduced due to the non-orthogonal shape of 
the grid regions that effect the transition from the in- 
ner Cartesian grid to the outer spherical grid (see Fig. 
4 in Ott et al. 2004) and by noise at the part in ^^10^ 
level in the EOS table interpolation. Since the resulting 
turbules execute more than twenty overturns during the 

For future calculations, the outer boundary radius should prob- 
ably be positioned still further out. 



initial phase of convective instability, and this convective 
phase reaches a quasi-steady state, the initial conditions 
and the initial perturbations are completely lost in subse- 
quent evolution. The seeds for the later shock instability 
are the non-linear convective structures that arise in the 
first post-bounce tens of milliseconds. Beyond these, we 
introduce no artificial numerical perturbations. 

3. SIMULATION RESULTS 

To demonstrate and represent the various core-collapse 
and post-bounce phases we think are important, we fo- 
cus in this paper exclusively on results using the 11- 
M© progenitor without rotation from Woosley & Weaver 
(1995), with the zoning and physics packages referred to 
in Progenitor dependences will be left to future work. 
The calculations were done from ~200 millseconds before 
bounce to ~660 milliseconds after bounce, significantly 
longer than most (though not all) multi-D simulations in 
the literature We have generated movies of the runs, 
which are available from the authors upon request. Here, 
we provide a sequence of stills depicting phenomena or 
transitions of relevance, along with analysis plots that 
help to clarify salient features. 

FiguresHthroughOlare an evolutionary sequence of the 
entropy distribution for the core region 750 kilometers 
on a side, with velocity vectors superposed to map out 
the evolving and complicated fiow fields. The sequence 
includes snapshots at times of 50, 150, 275, 310, 385, 
420, 470, 515, 575, and 608 milliseconds after bounce. In 
the last two stills, the supernova has clearly started to 
explode. 

In the sequences depicted in Figs. Hill we don't 
show the collapse phase because it is canonical and spher- 
ical. At ^^50 ms after bounce the shock has stalled, is 
roughly spherical, and is at a spherical radius (i?) of ^115 
km. Neutrino-driven convection has begun in the region 
~50 km wide interior to the shock wave. Despite the 
advection correction to the standard convective instabil- 
ity criterion pointed out by Fogiizzo, Scheck, & Janka 
(2005), we find that standard neutrino-driven convection 
at this stage has not been suppressed, though there is ev- 
idence of a weak £ — I pulsation similar to what has been 
identified by Fogiizzo and collaborators and by Blondin, 
Mezzacappa, & DeMarino (2003) (Figs. HandEJ. By 
^150 ms, the average shock radius has reached ^150 
km, and the convection is encompassing the region down 
to i? ~ 75 km. In the full angular region of 180°, we 
see 5-6 dominant turbules (eddies) with angular scales of 
~30-35°. Near ^200 ms, the average shock radius has 
receded back to ^110 km (Figs. 01and|5l). If shock re- 
cession had been our criterion for failure, we might have 
stopped the calculation here. However, at around ^^200 
ms, the shock is beginning to wobble up and down per- 
ceptibly in an ^ = 1 mode with a period near 25-30 ms 
and a AR/R near 25%. At ~250 ms, the AR/R is ap- 
proaching ^50% and the up-down asymmetry is quite 
pronounced. The growth time for the shock anisotropy 
varies, but is near 50-100 ms. We identify this early 
quasi-periodic oscillation with the advective- acoustic in- 
stability and the "standing accretion shock instability" 
(SASI) suggested by Fogiizzo (2001,2002), Foghzzo & 

^ Typically, one complete run will take four weeks on a state-of- 
the-art Beowulf cluster. 
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Tagger (2000), Foglizzo, Galletti, & RufFert (2005), and 
by Blondin, Mezzacappa, & DeMarino (2003), referred 
to from now on as the "shock instabihty." 

By ~300-350 ms, AR/R has grown to a factor of two 
and the wobble is taking on a more vigorous character. 
However, by ~275 ms it has clearly entered the nonlinear 
regime. In particular, when one side is executing its out- 
ward oscillation, the material in its outer extent reverses 
its flow direction out of phase with the material further 
in. The result is the generation, where these two regions 
collide, of a secondary shock wave; the flow now has 
nested shock waves. The creation of secondary shocks 
serves to saturate the amplitude of the shock oscillation, 
which nevertheless continues. Furthermore, the corruga- 
tion of the outer shock, where there is a kink in the shock 
normal, propagates in latitude in both hemispheres and 
collides at the poles, first at 0° then at 180°. In part, this 
shock focussing is an artifact of the 2D nature of the sim- 
ulation, but similar effects are seen in 3D pure hydro sim- 
ulations (J. Blondin, private communication) and are ex- 
pected from generic considerations. Each time the shock 
kinks collide, entropy and pressure {p) are generated and 
this pressure pulse moves at the speed of sound to the op- 
posite side. Importantly, at the shock intersection kinks, 
clearly seen in Figs, ^andl^lat 310 and 385 ms, the accre- 
tion through the outer shock is channeled and penetrates 
in lower entropy streams into the interior and onto the 
core. This phenomenon was previously seen by, among 
others, Herant et al. (1994), Burrows, Hayes, & Fryxell 
(1995), Fryer & Warren (2002), Rampp & Janka (2002), 
Scheck et al. (2004), Janka et al. (2005ab), and Buras 
et al. (2005). Though AR/R can reach factors of 3 or 
greater, this ratio always eventually decreases, if tem- 
porarily, and stabilizes. The angle-averaged outer shock 
radius does not increase monotonically. For instance, 
though at ~385 ms, one lobe of the outer shock reaches 
-400 km, as Fig. El indicates, at -420 ms the average 
radius has receded temporarily to —250 km. The 2D hy- 
dro calculations of Blondin, Mezzacappa, & DeMarino 
(2003) that identified this instability were done without 
neutrino physics, and with an artifical inner boundary; 
they were done without electron capture, neutrino heat- 
ing, neutrino cooling, and the core. Nevertheless, with 
neutrino physics and the core fully included, we verify 
the existence of this shock instability, first promoted by 
Foglizzo (2000,2001). 

However, we have discovered significant and important 
differences with the work of Blondin, Mezzacappa, & De- 
Marino (2003) and with the expectations from the work 
of Foglizzo and collaborators. First, the shock instabil- 
ity does not lead to an increase in the average shock 
radius in the first 150-200 ms after bounce. Due to elec- 
tron capture at the shock and to the quasi-hydrostatic 
sinking of the inner core caused by steady neutroniza- 
tion and cooling, an = 0" component is not in evi- 
dence. Shock asphericities with some £ = 1 component 
are generated, but mostly by the turbulent motions of 
neutrino-driven convection during this early phase. It is 
only after the core has settled, the electron capture rates 
at the shock have decreased due to the gradual decrease 
in the density of the accreted material, and there is no 
net energy loss behind the shock in the gain region that 
the advective-acoustic behavior identified by Foglizzo is 
clearly manifest. This delay could not have been cap- 



tured in the calculations of Blondin et al. (2003). Sec- 
ond, when the shock instability can finally be identified, 
the matter is already convecting nonlinearly due to neu- 
trino heating. Hence, a linear growth analysis may be 
inaccurate, or, at the very least, the seed perturbations 
for the shock-instability are the turbules and plumes of 
neutrino-driven convection. Third, the nonlinear phase 
of the shock instability brings with it secondary shock 
waves and the shock oscillation becomes saturated and 
non-periodic. By itself, the shock instability is not lead- 
ing to explosion and the average radius of the outer shock 
ceases to increase, though the wobble and top-bottom 
asymmetry can still be extreme. Fourth, within the first 
100-150 ms of the shock instability its behavior not only 
ceases to be periodic with a clear oscillation period, but 
is also not a simple normal mode. A new phenomenon 
arises. The increasing vigor of the turbulence in the inte- 
rior outside of the core (< 100 km), in part generated by 
the fluctuations and "wagging" over the core's surface of 
the lower-entropy accretion streams stirring this inner re- 
gion, begins to generate traveling acoustic waves, sound 
that propagates outward, with initial periods <5— 10 ms 
determined by the turbulc overturn times in the interior 
and the timescale for the heads of the accretion plumes 
to traverse the latitudes of the PNS surface. Such acous- 
tic power is a new feature in supernova theory. This is 
not the acoustic component of the "advective-acoustic" 
instability, which is a normal mode with a period near 
the sound travel time across the shocked region of —25- 
30 milliseconds. Rather, it is the propagating acoustic 
flux whose source begins as the turbulent energy in the 
roiling interior, driven by aspherical accretion of mat- 
ter and entropy onto the core. This is analogous to the 
"forced-turbulence" source for the acoustic flux studied 
by Goldreich & Kumar (1988,1990) in the context of 
solar convection. These authors showed that the effi- 
ciency of sound generation in a driven turbulent region 
increases roughly as the Mach number {Al) cubed ^ . The 
Mach numbers seen in our simulation are on average 
steadily rising, evolving in the inner —100 km from —0.1 
to as high as —0.9. This turbulence is initially fed by 
the grossly aspherical accretion flows associated with the 
shock instability and the dancing accretion streams reg- 
ulated by it. 

However, when the acoustic flux becomes appreciable, 
it significantly modifies the shock instability. In our 
calculation, over a period from —350 to —450 ms, the 
motions of the shock and accretion flow transition from 
being dominated by the shock instability to being domi- 
nated by the core's acoustic flux. The propagating sound 
waves issuing from the core modify the outer flow to such 
an extent that the shock motion is no longer even quasi- 
periodic and the momentum and energy flux of the sound 
start to determine the flow pattern. In particular, the 
acoustic flux is refracted and reflected by the accretion 
streams and emerges very anisotropically, and more on 

^ Furthermore, they showed that the peak in the acoustic fre- 
quency spectrum is near H/v, which is close to what we witness at 
this stage in our simulation, where H is the characteristic turbule 
size and v is the characteristic turbulent speed. However, another 
characteristic driving timescale is the period over which the ac- 
cretion streams "dance" / "flap" over the surface of the PNS. This 
timescale is a bit longer than H/v until the core g-mode oscillation 
starts to dominate, at which point the flapping all but ceases. 
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one side of the core than another. At the same time, the 
accretion streams are pushed by the momentum of the 
acoustic flux more and more onto the other side of the 
core, whichever that happens to be. The positions of the 
accretion streams and the average direction of the acous- 
tic flux become anti-correlated. Figure O clearly shows 
the low-entropy accretion streams and their positioning, 
as well as the discontinuities in the velocity vector field 
that mark the emerging sound pulses. In fact, in prop- 
agating outward down the density gradient, the sound 
pulses steepen into shock waves, subsequent pulses catch- 
ing up with and reinforcing previous pulses. From ~400 
ms after bounce, sound waves make themselves felt more 
and more and the resultant nested shock waves become 
a feature of the early supernova phenomenon. 

By 500 ms after bounce, the acoustic power from the 
core is quite pronounced and is starting to power outflow 
and the beginnings of an explosion. Figure ^depicts the 
temporal evolution after 100 ms of the entropy profile 
along the poles. The position of the outer shock, the 
early decay of its radius, the wobble due to the shock 
instability after ^250 ms, the fact that the outer shock 
radii at 0° and 180° (the top and the bottom) are roughly 
out of phase until ~500 ms, the entropization due to suc- 
cessive pulses at late times, and the onset of explosion 
are all clearly seen. Figure depicts the evolution of the 
Mach number profile, also along the axis of symmetry 
and as a function of time. This figure recapitulates Fig. 
^ showing the shock oscillations and the onset of explo- 
sion, while also showing the growth of the Mach number 
with time and the secondary shocks at later times. Note 
that despite clear evidence on this figure for a periodicity 
in the Mach number of the accreted matter and, hence, in 
the flow from ^^100 ms to ~250 ms after bounce, the av- 
erage shock radius does not increase due to a shock insta- 
bility until afterwards. By ~500 ms, one lobe of the outer 
shock is near 1000 km and by '^660 ms it has reached 
^2300 km. The accretion streams are now perennially on 
one side of the core and the acoustic flux, while mostly 
emanating from the other side is radiating in all direc- 
tions, except up the accretion streams themselves. Fig- 
ureOdepicts this clearly, though at ^^600 ms some matter 
is still infalling, particularly the lower-entropy streams. 
The entropy of some of the shocked matter can now be 
hundreds of units, due mostly to progressive shock heat- 
ing by the steepening sound pulses, not the neutrinos. 
Starting around '^550 ms, as the explosion commences, 
entropization due to these multiple shocks is clearly seen 
in Fig. 21 

After ~500 milliseconds, and certainly by 550 ms, the 
sonic power, though its early source was the turbulence 
around the core region, is driven mostly by a core os- 
cillation that is being excited by the violent accretion 
streams. This oscillation can be seen as early as ^^350 ms 
after bounce and has a period near ~3 milliseconds, very 
much smaller than the sound-travel-time in the shocked 
cavity. It is predominantly an = 1 g-mode of the inner 
^40 kilometers that has grown strong over a period of 
^100 ms to reach nonlinear amplitudes by ~500 ms. This 
mode would have been suppressed had we excised the 
inner core, not performed the calculations over the full 
180°, or performed the simulations in ID interior to some 
convenient radius and is discussed in ^ The core oscil- 
lation is driven by the energy in the accretion streams 



and by the turbulence around the core, both of which 
ultimately derive their energy from the gravitational en- 
ergy of infall. The oscillation is damped by sound waves 
that emerge out of the core. These sound waves steepen 
into shock waves, and, by dint of their momentum and 
energy flux , "ignite" the supernova explosion. The 
core oscillation is acting like a transducer for the conver- 
sion of the gravitational energy of infall into radiating 
acoustic power and at the later stages is a far more im- 
portant source of acoustic power than the inner turbu- 
lence. Moreover, from ~400 ms to ~660 ms the efficiency 
for the conversion of accretion power into sound power 
is increasing. As long as the accretion continues during 
this phase, the core oscillation seems to be driven and the 
sound is emitted. After the explosion has progressed suf- 
ficiently and accretion subsides, the core oscillation de- 
cays and the sound source abates. It seems that as long 
as the acoustic power due to core oscillation is needed to 
drive the explosion, it is available. If the neutrino mech- 
anism does not abort this scenario by inaugurating an 
earlier explosion, this may be a natural self-regulating 
mechanism for the supernova phenomenon. 

4. ACOUSTIC WAVES AS THE CRITICAL POWER AND 
MOMENTUM SOURCE 

Hydrodynamic, neutrino, convective, viscous (Thomp- 
son, Quataert, & Burrows 2005), and magnetic (Akiyama 
et al. 2003) mechanisms for driving core-collapse super- 
nova explosions have all been proposed and investigated. 
Acoustic power is (or would be) a new paradigm. If the 
neutrino mechanism obtains, it needs to inaugurate ex- 
plosion early, perhaps in the first 200-400 milliseconds 
after bounce. This is because the driving neutrino lumi- 
nosity, the absorbing mantle mass, and the optical depth 
of this mass to the emerging neutrinos, are decreasing. 
The result is that the neutrino energy deposited in the 
gain region is inexorably diminishing. In the simulations 
we present here, we do not see a neutrino-driven explo- 
sion. For all the most detailed simulations by the Garch- 
ing group (cf. Buras et al. 2005), except for their model 
sl5Gio_32.a, which they exploded artificially (see also, 
Janka, Buras, & Rampp 2003), this is also the case. This 
does not mean that the neutrino mechanism does not ob- 
tain. Indeed, 1) the LANL group (cf. Fryer & Warren 
2002,2004) does get neutrino-driven explosions using its 
gray SPH code, 2) many calculations to date may have 
critical flaws, and 3) accretion-induced collapse may well 
explode by the neutrino mechanism rather easily. 

However, there are certain virtues to acoustic driving 
that bear mentioning. First is that while the acoustic 
luminosity is much smaller than the neutrino luminosity, 
almost all of the sound is absorbed in the mantle mat- 
ter. At late times in our simulation, less than a percent 
of the Ve and De neutrino luminosity is absorbed. This 
amounts to an neutrino absorption power of <10^'' erg 
s~^, compared with an estimated core acoustic power at 
the end of our calculation near ~10^^ erg s-i. Fi gureEl 
compares the i^e and De neutrino luminosities, the inte- 
grated neutrino heating rate (power) in the gain region, 
the gravitational accretion power {Eacc ^^^M), and 
an estimate of the acoustic power due to core g-mode 
oscillations. Due in part to the very complicated flow 

^" the latter, minus that fraction radiated away by neutrinos 
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patterns and sonic refraction and reflection, it is difficult 
to determine precisely this total acoustic power. How- 
ever, with a simulation- motivated estimate of ~120 ms 
for the acoustic damping e- folding time (te) of the £ = 1 
core g-mode and our calculated energy (Eg) in this 
normal mode as a function of time, we can estimate this 
late-time acoustic power (~ Eg/rE) and have included 
this estimate in Fig. |B1 There is a time, in our calcu- 
lation near ~350— 400 ms, when the absorbed neutrino 
power goes below the rising acoustic power. 

Second, sound carries not only energy, but momen- 
tum, and this factor seems to be important in our sim- 
ulations. The momentum flux for sound with the same 
energy flux as neutrino radiation is larger by the ratio 
of the speed of light to the speed of sound, which in 
the inner mantle regions is as much as a factor of ten. 
Third, acoustic power propagates from where it is gen- 
erated to where it is needed; it fulflUs the central re- 
quirement of a core-collapse supernova mechanism that 
it involve energy transfer from the bound interior PNS 
to the outer exploding mantle. If the acoustic power 
is large enough, it is the ideal transfer agent. Fourth, 
the acoustic source seems to grow just when the neu- 
trino luminosity is ebbing and, importantly, it continues 
until explosion ensues. Fifth, the successive merger of 
trains of sound waves that steepen into shocks provides 
a non-neutrino way to entropize some of the matter and 
naturally achieve r-process conditions. 

5. THE EXCITATION AND DAMPING OF CORE 
OSCILLATIONS 

Within ^200-300 ms of bounce, the inner core has 
reached a total baryon mass above ~1.3 M0. Exterior 
to a radius of ~40 km, the density profile is falling off 
precipitously. The entropy in this inner core is of or- 
der unity (in units of Boltzmann's constant per baryon), 
and the outer Ye profile has a steep negative gradient. 
The neutrinospheres reside in this low Ye regime at radii 
of ~30— 40 km. Surrounding this inner core is a high- 
entropy accreting, turbulent mantle, whose entropy and 
Mach numbers are increasing with time. At ^--^200 ms, 
the total accretion rate is ~0.1 Mq s""'^ and is decreasing, 
but the temporal and spatial variations in the accretion 
fluxes and the ram pressures of the infalling plumes that 
eventually settle onto the inner core are quite large. Fig- 
ure [7| shows the time evolution of the spherical harmonic 
weighting coefficients, a^, of a decomposition into Y^^s 
of the fractional fluctuations of the pressure field at a 
radius of 35 km for £ = 0, 1, 2, and 3^^. We see in Fig. [7| 
the emergence of a pronounced i = \ component, whose 
modulation up to ~600 ms is at the period of the outer 
shock instability (^^25— 30 ms). 

This situation is perfect for the excitation of the nor- 
mal modes of oscillation of the inner core that are anal- 
ogous to classical stellar pulsations (Goldreich & Kel- 
ley 1977). Stars can execute p-modes, g-modes, and/or 
f-modes. The dominant restoring force for p-modes is 
pressure and that for g-modes is gravity. The f-modes 
generally have no radial nodes and are fundamentally of 

One can use the FWHM of the power spectrum of the core 
pressure fluctuations around the 3— ms mode or the early growth 
rate of the core oscillation kinetic energy to get a handle on this 
quantity. 

12 Ap(e) =<p>Y.iaiym{9) 



mixed type, though all realistic pulsational modes have 
some mixed character. Importantly, all the prominent 
modes excited in our simulations are significant admix- 
tures of both g-type and p-type character. We have per- 
formed an analytic modal analysis of the PNS core for 
its structure at ^^500 ms after bounce and display the 
derived periods in Fig. |S1 The p-mode (green dots) and 
g-mode (red stars) branches for each spherical harmonic 
angular eigenfunction (Y^^) are shown here as a func- 
tion of . Given an for predominantly g-modes the 
period increases with number of radial nodes and for pre- 
dominantly p-modes the period decreases with number of 
radial nodes, the latter as expected for sound waves in a 
cavity. Generally, the shorter period (higher frequency) 
oscillation modes are predominantly p-modes and those 
at longer periods (lower frequency) are predominantly g- 
modes and for all modal types the period is a decreasing 
function of The radial eigenfunctions of p-modes peak 
on the outside, while those of g-modes peak on the inside. 
At the boundary between the p-modes and the g-modes 
on Fig. 121 the modes are quite mixed, so that some of 
those modes might be assigned differently by others, for 
instance as f-modes. 

The periods for most of the modes depicted in Fig. 
|S1 are shorter than the characteristic times of pressure 
fluctuations and turbule overturn in the turbulent outer 
regions. Moreover, the dynamical times in the core- 
bounding region from 50 to 100 kilometers are longer 
than those in the core itself. In addition, though exci- 
tation by a neutrino K-mechanism has yet to be inves- 
tigated, one can show that the high-frequency p-modes 
would damp so quickly by sound emission that they can 
not be excited to signiflcant amplitude. However, as Gol- 
dreich & Kumar (1990) have shown, the efficiency for tur- 
bulence to excite g-modes is linear in Mach number, M . 
As M grows, the fraction of the turbulent power around 
the core that is pumped into g-modes increases. Further- 
more, the g-modes have the longest periods among the 
modes that could possibly be excited, and are, hence, 
the most likely to be resonantly excited by convective 
turbules and the waving accretion streams. 

Therefore, we would expect a whole spectrum of core 
modes to be excited, with the g-modes eventually dom- 
inating. A variety of is above should be represented, 
but the fact that the accretion is channeled into a small 
number of plumes/streams (generally, in our 2D calcu- 
lations, one per hemisphere at late times) means that 
lower values of £ should prevail. This is what we see in 
our 2D simulations. After ~450 ms, over a period of 
~100 ms, the € = 1 dipole mode with predominantly g- 
type character at a period near ^^3 milliseconds emerges 
to signiflcant amplitude, reaching Lagrangian displace- 
ments of ^^3 km. This core mode has a radial node near 
6—10 km, though it has a few nodes of p-type character 
further out. This mode corresponds to the analytic mode 
circled in Fig. It damps by the emission of acoustic 
power, sound waves, with an oscillator Q value (= lote, 
where te is the e-folding time for energy loss) near ~200. 

Without rotation, we expect the modes to be degenerate in 

m. 

l'* There are no ^ = g-modes, and no g-modes without nodes, 
since this would not conserve momentum. Note that the £ = I p- 
mode without a radial node would be a zero- frequency core trans- 
lation, a "kick." 
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In its nonlinear phase after ^580 ms, the total pulsation 
energy is '^10^'^ ergs and the radiated acoustic power is 
near ~10^^ erg s~^. Figure IHI depicts the time evolution 
of the spectrum of pressure fluctuations at a radius of 
30 km. On this frequency-time plot one sees the emer- 
gence of the 3— ms core g-mode oscillation to dominate 
after '^450 ms, as well as a number of other periodicities 
(modes) of lesser strength. For instance, the mode seen 
in Fig. I^lnear ~675 Hz is an = 2 harmonic. 

Though the radiation pattern we actually observe 
is quite anisotropic and variable, and the sound is 
severely refracted and diverted by the complicated ac- 
cretion/turbulent flow patterns, to gain insight one can 
estimate the dipole acoustic power (E^^^) radiated by a 
sphere of radius R executing periodic linear translational 
motion (wobbling) of amplitude AR with period P in a 
gas of density p and sound speed c^. This is given by 
(Landau & Lifshitz 1959): 
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where pi2 — gcm^'^, k {— uj/cs) is the wave num- 

ber of the sound generated, uj is the angular frequency of 
the oscillation {— 2tt/P), and uq = ujAR. The numbers 
given in eq. ^ are representative quantities for the prob- 
lem at hand and the impedance term {{kR)* / {A+ {kR)"^)) 
has been set to 1/5. Since in reality there are complicated 
density and sound speed profiles, the application of eq. 
l|T]l is no more than a very crude approximation. Never- 
theless, we see from simple analytics that even for small 
values of AR the acoustic power radiated by core oscilla- 
tions can be quite large, and can be of importance in the 
supernova context. At high oscillation frequencies and 
surrounding pressures, even small amplitude oscillations 
generate competitive acoustic radiation. This is what we 
see in the more detailed simulations, in which we derive 
acoustic powers due to core oscillations of ~10^^ erg s~^. 
Almost all of this power is absorbed in the mantle and 
is potentially available to the explosion. This is to be 
compared (see Fig. (SJ with the accretion power at late 
times of: 
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Core oscillations can therefore convert accretion power 
into acoustic power with some efficiency. In our base- 
line simulation at the same epoch, L;^ (total) is '^5 x 10^^ 
erg s~^ and L^^/^^ is ~2 x 10^^ erg s"""^. In contrast, 
L^^/^^XTi, (where is the gain- region optical depth 

to Vel^e neutrino absorption) is <10^° erg s^^. These 
power/luminosity comparisons, also depicted in Fig. |3 
highlight the problem, and its potential solution. 

Importantly, what we see in our simulations is the self- 
excitation of the core g-mode; the inner core plus in- 
ner accretion region is a self-excited oscillator. The in- 
teraction between the sound pressure and the accretion 
streams (see, e.g.. Fig. that results in an anticorrela- 
tion at late times between the angular positions of the 
accretion streams and the average direction of the emit- 
ted sound, and in the emergence of the £ = 1 mode, are 



clear features of our simulations. However, the flow fields 
are quite complicated and variable. What we have yet 
to determine is whether the sound modulates the excit- 
ing accretion streams to create a feedback loop with gain 
and a time delay, the classic context for a feedback am- 
plifier, quartz crystal oscillator, acoustic feedback oscil- 
lator, klystron, or edgetonc oscillation. This possibility 
requires further exploration, but is intriguing. 

The core oscillations are acting like a transducer for the 
conversion of the gravitational energy of accretion into 
acoustic power. The oscillation continues as long as the 
accretion continues, ensuring that acoustic power persists 
until accretion stops. This may be one means by which 
a supernova, and its energy, are self-regulating, given a 
progenitor structure. The halt of accretion is one man- 
ifestation of explosion. In fact, the accretion streams, 
sequestered predominantly on one side of the inner core, 
pump up the core oscillation and store energy in it like 
a battery or capacitor. Even after accretion subsides the 
energy in the core motion will continue to be damped by 
acoustic radiation, neutrino viscosity, and gravitational 
2 j^^agliation, the former being dominant. However, as the 

gas ,pressure around the core decreases due to explosion, 

^ because of the increasing impedance mismatch between 
the core and the expanding mantle, the acoustic damp- 
ing time is bound to increase. We have yet to calculate 
this effect, but the core may oscillate for a "long time." 
However long it takes, most of the oscillation energy will 
eventually be "discharged" sonically into the supernova 
explosion. 

The angular distribution of the emitted sound is fun- 
damentally aspherical. The actual modal periods depend 
upon the nuclear EOS used, and will be different when 
general relativity is incorporated into the calculations. 
Relativity will decrease the period and increase the pres- 
sures that surround the core, but it will also decrease the 
core size. Therefore, how relativity will alter the acoustic 
power is not yet clear. Furthermore, how the accretion 
streams will behave in 3D, where the pattern of infalling 
plumes and its evolution are as yet unknown, remains 
one of the major uncertainties of our acoustic model. A 
rotation axis may set a direction and render important 
aspects of a 3D simulation like those in a 2D simulation. 
This is what is suggested in the linear analysis of Mi- 
ralles et al. (2004) and what one would conjecture from 
the generic alignment efi^ects on convective plumes of the 
S0lberg-Hoiland instability. Otherwise, there may be a 
spectrum of modes that are not degenerate in azimuth 
and in "magnetic quantum number" m. In addition, 
though we did not see much power in € = 2 modes, they 
may contribute more in the real 3D situation. Neverthe- 
less, from our calculations, uniquely allowing as they do 
core motion on a 2D grid, the excitation of non-radial 
core g-modes seems to be a natural consequence of the 
later stages of stellar collapse. 

6. THE ROLE OF NEUTRINOS 

In the baseline simulation we have performed for this 
paper, despite the fact that it involves an 11-Mq pro- 
genitor with a steep outer density profile and despite the 
appearance and strengthening of neutrino-driven convec- 
tion in the first ^200 ms after bounce, we do not witness 
a neutrino-driven explosion. This does not mean that 
the neutrino-driven explosion mechanism is ruled out. It 
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may be that different input physics, different progenitors, 
an improved numerical technique, or not employing the 
MGFLD approximation could result in a neutrino-driven 
explosion. Since they should accrete their outer bound- 
aries quickly within less than '^100 ms after bounce, 
ONeMg white dwarfs that achieve the Chandrasekhar 
mass and collapse due to accretion from a companion 
(accretion-induced collapse, AIC) can't avoid exploding 
in the early post-bounce stages before the shock instabil- 
ity and core acoustic sources turn on. Nevertheless, with 
our fully 2D multi-group code, general Poisson solver, 
and adequate energy-group and spatial resolution, shock 
stagnation is prolonged until long after the first peak in 
the average shock radius (near 150-200 km) and its sub- 
sequent initial decay. In the past, this peaking and decay 
were taken as signatures of a fizzle. 

However, if our result is an adequate representation 
of Nature, then what is the role of neutrinos in the 
more-delayed sound-driven mechanism we have discov- 
ered? First, neutrinos still deposit energy in a gain re- 
gion. Figure depicts not only the neutrino luminosi- 
ties, but the neutrino power deposited in the gain re- 
gion as a function of time. Hence, neutrinos will cer- 
tainly contribute a fraction of the explosion energy. Sec- 
ond, from the work of Foglizzo and collaborators and 
of Blondin, Mezzacappa, & DeMarino (2003), we expect 
that the shock instability requires a large shock height 
and a large ratio between the speeds of sound at the core 
and at the shock. Neutrino-driven convection results in 
larger shock radii after the stalling of the shock than ob- 
tain without neutrino-driven convection. Nevertheless, 
the shock instability can't commence until the inner core 
slows its quasi-static contraction, a phenomenon slaved 
to neutrino and lepton losses. Third, neutrino heating 
of the accreted and shocked material that goes through 
the gain region is one reason the entropies around the 
core increase after ^^200 ms after bounce. The other 
reason is enhanced, multiple shock heating due to the 
shock instability. Higher entropies and larger entropy 
perturbations result in larger overpressures during the 
earlier turbulence-driven acoustic flux phase seen in our 
simulations before the core oscillation becomes vigorous. 
Fourth, neutrinos are responsible for the failure of the 
direct mechanism for progenitors less massive than '-~^20 
M0^^. If this mechanism succeeded, then the gravita- 
tional mass of the residual neutron star very likely would 
be too small to explain observed pulsar masses. Fifth, 
the absorption of v^. and neutrinos by the ejecta may 
be responsible for ensuring that its Ye is not too small to 
be consistent with nucleosynthesic constraints (Frohlich 
et al. 2005; Pruet et al. 2005). Furthermore, YgS above 
0.5 that are possible only due to an excess of Ve absorp- 
tion over Vf, absorption might be needed to explain the 
rp-process isotopes. However, we find that the mass of 
the core ejecta in our acoustic-power-driven explosion is 
less than 0.01 and that there may not be as much 
of a problem with anomalously neutron-rich ejecta for 
the lowest-mass massive stars. Finally, neutrino cooling 
of the multiple shocks emanating from the core slightly 
diminishes the acoustic energy that in our scenario ulti- 
mately powers the supernova. Determining the correct 

For progenitors more massive, nuclear photodissociation alone 
can stall the shock. 



magnitude of this negative effect will be important and 
might require better neutrino transfer and better weak- 
interaction physics than we have incorporated into our 
2D MGFLD approach. 

In the calculations we present in this paper, neutrino 
transport and weak-interaction physics are crucially im- 
portant, but arc not the critical factors in the explosion 
mechanism. Anisotropic neutrino radiation pressure may 
in some circumstances 1) affect the core oscillations, and 
2) contribute an integrated net impulse to the PNS and, 
therefore, a kick. However, in our baseline model, with 
the single progenitor on which we focus, such effects are 
subdominant and the direct role of neutrinos is subordi- 
nate to the other factors described in inland SjSl 

7. DISCUSSION AND CONCLUSIONS 

In this paper, we present a new mechanism for core- 
collapse supernova explosions that focusses on the acous- 
tic power generated in the core region. The strength, 
radiation pattern, and character of the emergent sound 
are influenced by the shock instability that arises after 
^200 ms, but this instability is not the agent of explo- 
sion. Rather, it is the acoustic power generated first 
in the turbulent region around the inner PNS core and 
then through the excitation and sonic damping of core 
g-mode oscillations. An £ = \ mode grows at late times 
to be prominent around ~500 ms after bounce, though 
it is in evidence as early as ~300 ms after bounce. In 
our Newtonian simulation, its period is ~3 ms. At the 
end of the calculation, this core g-mode contains '^10^° 
ergs, is radiating sound into the exploding mantle at a 
rate near ~10^"'^ erg s~^, and has a Q- value of ^200. At 
late times, but before explosion, the accreting protoneu- 
tron star is a self-excited oscillator, "tuned" to the most 
easily excited core g-mode. After the acoustic power be- 
comes strong, the average direction of the radiated sound 
and the angular positions of the exciting plumes out- 
side the core are anti-correlated. The possibility that 
the PNS-core/accretion-stream system is a feedback am- 
plifier is an intriguing notion suggested by our simula- 
tions, but that has yet to be properly explored. Since 
core H. = Q modes are not excited to an appreciable de- 
gree, the driving acoustic radiation pattern is fundamen- 
tally anisotropic and the initial phase of the explosion 
is unipolar. In addition, due to the extreme breaking 
of spherical symmetry, this model manifests simultane- 
ous accretion and explosion. The ultimate source of the 
acoustic power is the gravitational energy of infall and 
the core oscillation acts like a transducer to convert this 
accretion energy into sound. 

An advantage of the acoustic mechanism is that acous- 
tic power does not abate until accretion subsides, so that 
it is available as long as it is needed. This is not the 
case for neutrinos, whose luminosities and mantle heat- 
ing rates inexorably decrease at late times. Hence, this 
may be the long-sought-after self-regulating mechanism 
of the supernova energy and, being a function of the ac- 
cretion rate, is determined mainly by the progenitor den- 
sity structure, with some ambiguity due to the sensitivity 
to initial conditions in chaotic flows and due to rotation. 
Furthermore, unlike neutrinos, sound pulses that steepen 
into shock waves are almost completely absorbed in the 
baryonic outer mantle. Moreover, sound deposits not 
only energy, but momentum. Ironically, the large accre- 
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tion rates that were thought to inhibit explosion in the 
neutrino-driven mechanism are not a disadvantage in the 
acoustic model. In fact, the very accretion that might be 
inhibiting the neutrino mechanism facilitates the acoustic 
mechanism, but at later times. Furthermore, by their na- 
ture the accretion funnels constitute a disproportionate 
share of the infalling material. The result is that while 
at the same time these streams are exciting the core os- 
cillations accretion from the other directions is weaker, 
thereby presenting less of an obstacle in those directions 
to eventual explosion. In our calculation with the 11- 
Mq progenitor, the delay to explosion is '^550 ms. This 
is longer than the delay traditionally associated with the 
onset of multi-D neutrino-driven explosions (^200— 300 
ms). 

The baryon mass of our remnant neutron star is ^^1.42 
Mq, with a gravitational mass near ~1.3 Mq. This is 
close to what is canonically expected from measured val- 
ues. Had the core exploded much earlier, the mass re- 
maining might have been uncomfortably lower. 

The successive superposition of multiple shocks and 
the shock-shock interactions during the shock instabil- 
ity phase lead naturally to higher entropies in a fraction 
of the supernova ejecta. In our baseline calculation, en- 
tropies higher than 300 were achieved. High entropies 
may be a partial consequence of the thin accreting enve- 
lope, suggesting that higher-mass progenitors would not 
so easily yield such high numbers. Since the supernova 
energy, kick velocities, and ejecta entropies are all depen- 
dent in our acoustic model on the accretion regime, these 
quantities may in fact be correlated along the progenitor 
mass continuum, with lower-mass progenitors having low 
explosion energies, some r-process nucleosynthesis, and 
weak kicks, and higher-mass progenitors having higher 
explosion energies, no r-process, and significant kicks. 
Observations do suggest that lower mass massive stars 
are the likely sites of the r-process (Mathews, Bazan, & 
Cowan 1992) and that supernova energies might span a 
wide range (Hamuy 2003). Be that as it may, the progen- 
itor dependence of the supernova explosion systematics 
is an important aspect of the supernova story and one 
that we have yet to explore in the context of this new 
mechanism. 

Though we find that an ^ = 1 g-mode eventually dom- 
inates, ^ = 2 and ^! = 3 modes and harmonics are in evi- 
dence and there is also likely to be nonlinear mode-mode 
coupling. Unlike i = \ modes, 1 — 2 modes will generate 
gravitational radiation and will do so at characteristic 
frequencies (!) that is a function of the EOS, relativ- 
ity, and the PNS structure. Using the results from ear- 
lier quasi-static PNS cooling calculations, Ferrari et al. 
(2003) calculated its g-mode frequencies and their evo- 
lution. Including as they do general relativity, Ferrari et 
al. derive frequencies for the fundamental g-modes that 
are slightly higher than we derive using our Newtonian 
code, and find that these frequencies, after an initial rise, 
evolve to lower values. Our calculations recapitulate ex- 
actly this initial rise (Fig. El and subsequent decay after 
~670 ms (not shown in Fig. The total energy ra- 
diated in gravitational waves from this non-rotating 11- 
Mq model is ~10~^ Mq equivalent (Ott et al. 2005, in 
preparation). Hence, in the excitation of normal modes 
in the supernova context, we may have a direct signature 
of core physics and supernova phenomenology. 



Our calculations are Newtonian and were performed in 
2D. Relativity and a different nuclear EOS (Ferrari et al. 
2003) will change the frequencies of the core modes and 
the efficiency of sound generation, the latter in ways that 
are not yet obvious. Three-dimensional simulations are 
likely to be different from 2D simulations, manifesting 
more realistic plume structures. The 2D/3D difference 
is the major unknown and uncertainty in our scenario. 
However, it would be difficult to suppress in 3D the ex- 
citation of core g-modes and the generation of acoustic 
radiation that in 2D becomes so pronounced. Further- 
more, even slight rotation may set a natural axis for the 
evolution in 3D. Otherwise, more subtle initial conditions 
will break the symmetry. 

We have also performed a simulation that restricts the 
hydrodynamics to a 90-degree wedge, with a reflecting 
boundary at the equator, as opposed to operating in 
the full 180-degree domain. In this case, we see that 
the £ = 1, 25-30 ms wobble due to the shock instabil- 
ity is completely suppressed. However, instead of aris- 
ing near 200-250 ms after bounce, the shock instability 
starts later, around 450-500 ms after bounce, and in an 
t — 2 mode. Before that, the shock maintains a roughly 
spherical shape. Though the delay is much longer, the 
shock executes small-amplitude, long-period £ = 2 mo- 
tions until roughly ~650 ms after bounce, at which time 
the shock instability becomes vigorous. Just before that, 
near ^600 ms, the PNS starts to radiate sound waves 
from a predominantly i = 2 core g-mode with a ^400 
Hz frequency (see Fig. IHJ. Clearly, 90-degree calcula- 
tions suppress an important component of the supernova 
story, but just as clearly, the same instability and oscil- 
lation phenomena, though dominated hy l = 2 morphol- 
ogy, eventually arise. It just takes longer for 1 = 2 modes 
to grow. 

It is important to list the reasons the acoustic phenom- 
ena we have identified and presented in this paper were 
not seen before. First, most calculations were stopped 
after the shock radius first subsided around 200-300 ms 
after bounce, but before the shock instability was much 
in evidence, and before turbulence around the core could 
generate significant acoustic power. Second, those calcu- 
lations that were not stopped early were continued be- 
cause they experienced an early neutrino-driven, multi-D 
explosion. Such an explosion arose either naturally from 
the particular code being used, or was artificially pro- 
duced. If the explosion commences early, the PNS core 
oscillations are not excited to useful amplitudes and the 
shock instability is more mild. Third, and most impor- 
tantly, to date all other grid-based supernova codes have 
conducted calculations either with the cores excised, han- 
dled in ID, or on a ~90° wedge, thereby completely sup- 
pressing core oscillations and the resulting ^ = 1 acoustic 
flux. Such procedures can even muffle the acoustic flux 
generated in the turbulent inner "convective" zones. As 
a result of some combination of the reasons above, no 
previous supernova simulations, before those using VUL- 
CAN/2D, could have discovered the acoustic driving and 
core oscillation mechanism. 

Therefore, one key to the discovery of this poten- 
tially important mechanism was the computational lib- 
eration of the inner core to execute its natural multi- 
dimensional motions. Another key was patience to per- 
form the simulations to very late times. It may be that 
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better neutrino and weak-interaction physics, full 3D ra- 
diation/hydrodynamic simulations, multi-group/multi- 
angle simulations, a new suite of progenitor models, the 
use of other massive-star progenitors (we have here stud- 
ied only one), or some qualitatively important flaw in 
our approach or implementation will alter our conclu- 
sions here. The role of rotation must be explored, for 
both the neutrino and the acoustic models. Interestingly, 
when magnetic fields are incorporated into the calcula- 
tions, the generation of Alfven waves at the core oscilla- 
tion frequency might prove to be another power source 
for the supernova (S. Woosley, private communication). 

The multi-D neutrino-driven mechanism may still ob- 
tain, and likely does so for AIC systems. Moreover, it 
may be that some supernovae explode by the neutrino 
mechanism^^, while others, if the neutrino mechanism 
fizzles, explode by the acoustic mechanism or an MHD- 
jet mechanism (Akiyama et al. 2003), though the latter 
requires very rapid rotation that may not be available in 
the generic core-collapse supernova context"'^^. However, 
though unlikely to be correct in detail, our calculations 
and analysis may be pointing to important new phenom- 
ena in the theory of core-collapse supernovae that have 
hitherto been overlooked. Much remains to be done, in- 
cluding the determination of the progenitor-mass and 
-profile dependence, explorations in 3D with neutrino 

^® See the papers by Janka et al. 2005ab for the suggestion 
that the neutrino mechanism might still obtain for the lowest mass 
massive stars when the full 180° computational domain (even with- 
out liberating the core to execute oscillations) and state-of-the-art 
Boltzmann transport are incorporated. These authors seem to 
see the onset of explosion for an 11.2-Mq model from Woosley, 
Heger, & Weaver (2002). Interestingly, this model has a small iron 



transport, and verification and scrutiny by other groups. 
Nevertheless, we believe that the new lines of investiga- 
tion we have opened up in this paper, and the potential 
of the acoustic mechanism interpreted broadly, should 
stimulate many future studies in this central problem in 
theoretical astrophysics. 



We thank Rolf Walder, Thierry Foglizzo, Todd Thomp- 
son, Itamar Lichtenstadt, Stan Woosley, Casey Meakin, 
Chris Fryer, and John Blondin for fruitful discussions 
and their insight. We acknowledge support for this work 
from the Scientific Discovery through Advanced Comput- 
ing (SciDAC) program of the DOE, grant number DE- 
FC02-01ER41184, and from the NSF under grant AST- 
0504947. E.L. thanks the Israel Science Foundation for 
support under grant # 805/04, and C.D.O. thanks the 
Albert-Einstein-Institut for providing CPU time on their 
Peyote Linux cluster. This research used resources of the 
National Energy Research Scientific Computing Center, 
which is supported by the Office of Science of the U.S. 
Department of Energy under Contract No. DE-AC03- 
76SF00098. Finally, we thank Don Fisher, Youssif Al- 
nashif, and Moath Jarrah for their help generating both 
the color stills and the movies associated with this work. 

core mass of '^1.25 Mq, a density cliff near 1600 km, and a much 
steeper mantle density gradient around 1000-2000 km than the 11- 
Mq model from Woosley & Weaver (1995) upon which we focus in 
this paper. 

However, such rapid rotation might enable the MHD scenario 
for hypernovae and long-duration gamma-ray bursts. 
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Fig. 1. — Colomiap stills of the entropy, taken at 50 (top left), 150 (top right), 275 (bottom left), and 310 ms (bottom right) past core 
bounce, with velocity vectors overplotted. Here "Width" refers to the diameter through the middle; the radius through the middle is 375 
kilometers. Note that on this figure, as well as on Fig. |^ for ease of comparison between panels the same colormap is used. It extends 
up to 100 units (red), above which it saturates (See text for discussion). These calculations have been done for a full 180° and the axis of 
symmetry is vertical. 
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Fig. 2. — Same as Fig.^ but this time showing the entropy at 385 (top left), 420 (top right), 470 (bottom left), and 515 ms (bottom 
right) past core bounce (See text for discussion). 
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Fig. 3. — Same as Fig. but this time showing the entropy at 575 (left) and 608 ms (right) past core bounce. Note the acoustic waves 
emanating from the core, most easily seen in the velocity vector field. The color map extends to entropies of 200 (red), and then saturates 
for entropies beyond 200. The low-entropy accretion streams that are exciting the core g-mode are clearly seen (see text for discussion). 
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Fig. 4. — Time evolution, from 100 ms to ~660ms after core bounce, of the entropy (logarithmic scale) along the axis of symmetry, i.e., 
at r = or latitudes ±90°, and covering the inner 2500 kilometers of the grid. An entropy ceiling of 100 has been adopted. 



16 




0,1 D.2 0.3 0.4 0,& 0,6 

T\mB after bounce {s) 



Fig. 5. — Time evolution, from 100 ms to ~660 ms after bounce, of the Mach number along the axis of symmetry, i.e., at r = or latitudes 
±90°, and covering the inner 2500 kilometers of the grid. Note how clearly the multiple secondary shocks emerging at late times can be 
discerned and that the displayed Mach number saturates at a value of two. 
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Fig. 6. — Time evolution, after core bounce, of the sum of the v^^ and neutrino luminosities (Live + t'e), red line); the gravitational 
energy of the accreted material (blue) through a radius R =35 km, defined as (dE/dt)i^cc = ^^^Mn, where Mji is the mass interior to 
R and Mn the mass accretion rate through R; the net energy deposited by neutrinos in the gain region (green), defined as {dE/dt)^^i^ = 
Jf^ e(R, 9)dm{R, 8), where Q is the gain region, e{R, 9) is the net neutrino heating rate, and m,{R, 9) is the mass of the cell at (i?, 6); and 
the acoustic power ((cfS/dt) acoustic) radiated by the core oscillation. This is approximated here by the ratio Eg/r^, where Eg is the total 
g-mode energy interior to a radius of 40 km and te is the decay time of the core pulsation, taken as 120 ms. The default acoustic power 
(solid black line) threads through a ±50% swath that represents our estimate of the current ambiguity in extracting the numerical acoustic 
power from the simulation. 
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Fig. 7. — Time evolution of the spherical-harmonic coefficients for the fractional pressure variation for modes I =0 (black), 1 (red), 2 
(blue), and 3 (green) at a radius R = 35 km, given by at = 2tt de sineY^{e){P{R, 6)- < P{R,e) >g)/ < P(R,e) >g. Notice that 
despite the fact that the £ =1 mode looms large, the £ =2 and £ =3 modes are also in evidence. The £ =2 (harmonic) mode will result in 
a distinctive signature in gravitational radiation detectors, initially at a frequency near ~675 Hz. This frequency is likely to be different 
(higher) when general relativity is included (Ferrari ot al. 2003). 
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Fig. 8. — Frequency, v, versus the index, I, of the spherical-harmonic functions o(^i<^) foi' the analytic core g-modes (red stars) and 
p-modes (green circles) obtained from a spherical average of the full 2D simulation profiles at 500 ms after bounce. The corresponding 
periods (P = are given on the right axis. Nearly all modes have some g-type nodes and p-type nodes. The £ = 1 mode highlighted 

with the circle corresponds to the mode with predominantly g-mode character which has been most easily excited in our simulation. 
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Fig. 9. — Colorscale of the power spectrum of the fractional pressure variation {P{R, 8)— < P{R,6) >g)/ < P{R,6) >g at a radius 
R = 30 km, as a function of time and frequency. For each time t, a power spectrum is calculated from a sample of time snapshots covering 
titSOms, at a resolution of 0.5 ms. Note the emergence of power in the ^SSO Hz (=3 ms) g-mode, as well as the strengthening £ =2 
harmonic mode near ~675 Hz at late times. The latter is of relevance for gravitational radiation emission. 



